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The recent measurements by CLEO for the inclusive J/*]/ and 5'(2S') production in T(IS') decay 
and our previous calculation are analyzed. The J/^' momentum spectrum and the production ratio 
of ^'(25') versus J/*!/ favour T J /^{^{2S)) + ccg as the dominant contribution. We point out that 
the differences between the experimental data and our previous results are mainly originated from the 
setting of the parameter charm quark mass. Encouraged by this result, we investigate the associate 
strange particle enhancement as a probe for the open charm particles in T — > J/'I'('I'(25')) + ccg. 
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Besides the important role in the electroweak and CP 
violation research, heavy quark (b or c) physics is also the 
traditional arena for Quantum Chromodynamics (QCD), 
both in perturbative (P) and non-perturbative (NP) as- 
pects. The production and annihilation of heavy quarks 
via strong interactions can be calculated by PQCD for 
their large masses naturally set the hard scales, while 
the structures of heavy hadrons are ruled by NPQCD. 
Since heavy quarks can hardly be produced in NPQCD 
phase, one can compare the quantities of them calculated 
by PQCD with those of the corresponding heavy hadrons 
to study the NPQCD effects. Hence, heavy quarks act as 
bridges linking the P and NP phases of QCD. Neverthe- 
less, such a distinct character requires a careful investi- 
gation on both PQCD and NPQCD phases for a certain 
process, as one has to do in the study of quarkonium 
production. 

Non-Relativistic QCD (NRQCD) [l| provides a factor- 
ization framework j^] to calculate the heavy quarkonium 
production and decay. For example, the inclusive produc- 
tion of charmonium He in the pp collision can be written 
as 



da{pp HiP)c+X) 



(1) 

which is a sum of various contributions. Each contri- 
bution is the product of two factors: One is the cross 
section of the quark pair production in a definite (colour, 
angular momentum, etc.) state, the other is the NRQCD 
matrix elements to describe the transition probability of 
the quark pair to the relevant quarkonium. This for- 
mula clearly implies that one has to take into account 
all the competing contributions in calculating the cross 
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section measured by experiments. For the inclusive char- 
monia production in pp interaction at Tevatron 3] , the 
colour-octet process [J| seems necessary to explain the 
total cross sections as well as the transverse momentum 
spectra (But it is not clear yet why the polarization of 
J/vP not properly described because only the cross 
section of all the colour-singlet processes summed up is of 
one order lower than the experimental data. This is one 
example that teaches us to include all kinds of contribu- 
tions, but does not mean colour-octet process(es) always 
the dominant. Until now, the values of colour-octet ma- 
trix elements can only be extracted via fitting the data 
but are generally believed to be smaller than the colour- 
singlet ones, based upon the NRQCD velocity counting 
rule. Furthermore, when the same matrix elements de- 
termined at CDF are applied to photoproduction of J/^ 
at HERA, the color-octet contribution is about a factor 
of ten too large This indicates that these matrix el- 
ements may be even smaller. So, for a concrete process, 
one can not judge a priori which kind of process is the 
main contribution and can not truncate the hard part (da 
in Equation (IT|)) to a certain order of as, without distinc- 
tion for colour-singlet or colour-octet partonic processes. 
The reason is that the suppression from as for higher 
order colour-singlet partonic process (es) can be compen- 
sated by the larger colour-singlet matrix element (s) (The 
colour-singlet matrix elements of vector mesons can be 
determined unambiguously from their leptonic decays). 
So higher order (in as) colour-singlet process could also 
contribute significantly. This is the key point to un- 
derstand the inclusive J/^ production data in T decay 
ii. 

The strong decay of T is dominantly via 3 gluons, be- 
cause of charge parity conservation. So the straightfor- 
ward partonic process for the charmonium production is 
T ggg* ggcc which is of lowest order PQCD 
to create charm quark pair. Here the charm quark pair 
is produced via the virtual giuon splitting and is inher- 
ently colour-octet. When the colour-octet matrix ele- 
ment extracted from Tevatron is applied, one can get the 
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JV* + ccg 
0. This is 



branching ratio comparable to the experimental data Q . 
However, because one of the real gluon is preferred to be 
soft, the J/4' momentum spectrum peaks at the maxi- 
mum value, like a 2-body final state configuration. This 
was not consistent with the old CLEO measurements [1] , 
but the error bar was too large to draw a final conclusion, 
then. Another colour-octet process which is also of the 
2-body final state configuration was suggested ig- 
norant of the J/^* momentum spectrum, too. Recently, 
the CLEO collaboration make a new measurement for 
the T(IS') decays to charmonium final state ^9]- The 
data sample is 35 times larger. The more precise mea- 
surement confirm a rather soft spectrum of the J/^* in 
the T J/'^ + X process, which is consistent with the 
colour-singlet process prediction ^ both in the branch- 
ing ratio as well as the J/'^ spectrum but clearly confiict 
with the colour-octet processes in J/^> momentum spec- 
trum. 

As is discussed in [Tj and in above, because the much 
larger colour-singlet matrix element (s) can compensate 
the suppression from powers of a^, one should not ig- 
nore the colour-singlet process (es) if the suppression from 
powers of as is not too large (e.g. here, only one ex- 
tra as). Hence we calculated the contribution of T — > 
ccg to the inclusive J/vt production of T decay 
is dominant of colour-singlet process and can be 
treated within the traditional non-relativistic wave func- 
tion approach [l^, [H, [l3l- For this process, it is very 
clear that the more final state particles, especially the 
open charm hadrons from the unbound charm quark pair 
make the J/'^ spectrum much more soft. That the domi- 
nant contribution of this channel is via colour-singlet pro- 
cess requires that the bound c and c should come from 
different virtual gluons and in nearly the same momen- 
tum. This configuration in phase space eliminates the 
2-body like configuration that the real gluon is soft and 
the 2 virtual gluon splitting into 2 back to back cc pair 
with minimum invariant mass. So the momentum spec- 
trum of J/^ can not peak at the largest value. Though 
the qualitative analysis is perfectly consistent with data, 
the new data [9| and our calculation [7| does not coincide 
completely. One may be interested whether there are 
other dynamical reasons ignored. Hence in the following, 
we first clarify that the "inconsistency" between the data 
and our calculation comes from the setting of the param- 
eter charm mass in numerical calculation. We emphasize 
on the J/vf momentum spectrum and the production ra- 
tio of ^{2S) versus J/^. Furthermore, we investigate 
the associate strange particle enhancement as a way to 
probe the charm quark accompanying with J/4'. 

The reason for the "harder" spectrum of our 
previous calculation than that of the experiment 

The CLEO collaboration claims that the momentum 
spectrum is "closer to, although softer than" our colour- 
singlet prediction 0. However, this is just because we did 
not take into account the open charm threshold effect in 
the partonic level calculation. In our calculation, we let 



the charm quark mass equals to half of J/4' mass. This 
leads to the unphysical region when the ccg cluster mass 
smaller than 2m d. If taking into account the threshold 
effect, the peak in our calculation should move to the left 
about 0.1 unit of x (scaled momentum). This means the 
position of the peak of our calculated J/^ momentum 
spectrum (see, FIG. 3 of or FIG. 10 of H) should move 
to around x = 0.4 rather than x = 0.5. Considering that 
the width of the bins of the experimental data is 0.2 unit, 
the prediction of the colour-singlet model is consistent 
rather than softer than the data. 

The "wrong" result of our calculation o/\E'(2S') 
production rate 

If we only consider the production rate of 4^ (25*) it- 
self, it may be natural to choose TOc = M^,(2s)/2 and 
use the 4'(25) e+e^ decay width to extract the wave 
function (In fact sometimes one needs to tune the quark 
mass to get the proper absolute value of the cross sec- 
tion). This is what we did in our previous calculation. In 
, when calculating the partonic cross section of charm 
pair creation for J/4' production, we use TOc = Afj/*/2; 
while calculating that for the 4' (25") production, we use 
TOc = -W*(2S)/2. As can be seen from the Table demon- 
strating the charm mass dependence, this leads to a dif- 
ference of about 2.5 times. The ratio of the wave func- 
tions square and the branching ratio of 4" (25") to J/4' 
lead to another factor about 4. This is why in our paper 
we only predicted about 10 per cent of J/4' are from the 
decay of 4* (25"). However, if we calculate the production 
ratio(es) between J/4', 4^ (25*), etc., and adopt that they 
are different eigenstates of cc bound state, their differ- 
ences should be only inherent in the wave functions. The 
M^^2S) is explained as originated from difference 

of binding energies. Hence the charm quark mass should 
be taken as the same. This makes the calculated results 
consistent and simple for the colour-singlet processes. In 
the framework of non-relativistic wave function approach 
[3 Ell 13 ! '^^'^ easily derive, for the S-wave parti- 
cles, 

B{T 4'(2S') + ccg) r(T ^ 4'(2S') + ccg) 



B{T J/4' -t- ccg) 



r(T ^ J/4' + ccg) 



r(4'(2g) -> e+e") 
r(J/4'^e+e-) ■ 



(2) 



In the above equations, '025(0) ^^"-^ V'is(O) are wave func- 
tions at origin for 4'(2S') and J/4', respectively. To get 
the relations in Equation ([2]), we notice that r(T 
4'(nS') -I- ccg) should be, to order 0(u°), 



l^,fs(0)l' 
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|M(T ^ c(P/2)c(F/2)cc5)|2, 



(3) 



where dR' represents other integral variables than the 
'^{nS) energy E for the '^{nS)ccg phase space ele- 
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ment. P is the four momentum of ^{nS). M{T 
c{P/2)c{P/2)ccg) is the invariant ampUtude for T 
c{P/2)c{P/2)ccg) process. In this formula, when taking 
the same charm quark mass, the only difference for dif- 
ferent S-wave bound states is the wave function square. 
For the Icptonic decay width of S-wave vector charmo- 
nia, a very similar formula and hence conclusion can be 
obtained. Only that one should take the energy of the 
charmonia at rest frame to be 2mc, neglecting their mass 
differences. 

Using the values of 2002 data book [l^, which is also 
used by the experimental group Q , one can calculate the 
ratio of the e+e~ decay widths (central value) between 
^'(25') and J/^-. It is striking to be 0.416 (the to-date 
values give 0.45), which is very near the central value of 
the data 0.41 However, the colour-octet process can 
not simply derive such a relation with the e+e^ decay 
width[20'] 

Unfortunately, we have not got the analytical formula 
for the partonic process J/'^ — > ccccg, so it is difhcult 
to evaluate the P wave particles which is also measured 
by CLEO. The calculation needs the derivatives of the 
partonic amplitude. Encouraged from the S-wave pro- 
ductions, we leave the calculation of P-wave particles as 
an interesting separate work to do. 

From the above discussions, it is clear not that 
NRQCD leads to the prediction that the colour-octet 
processes dominate in T J/^ + X. It is just the 
NRQCD formulae require a global consideration of all 
possible partonic processes to explain the data. Frankly, 
there are many uncertain parameters like quark mass, as 
(scale dependence is significant in the lowest order calcu- 
lation), etc. So only quantities like the relative produc- 
tion ratio and shape of spectrum which are not sensitive 
to the absolute values of these parameters are relevant to 
justify the partonic processes. Only to get a cross sec- 
tion or branching ratio comparable with data sometimes 
seems tricky. 

Associate strange particle enhancement The 
above arguments confirm the significance of the colour- 
singlet process. As has been pointed out by Q and the 
experimental group @, a very important probe for the 
partonic process is to measure the content of the cluster 
produced associated with the J/^ in the T decay, which 
is theoretical- framework- independent. It is interesting 
to notice that in e^e~ continuum process at nearly the 
same energy, BELLE Collaboration [l3| reported that 
about 60 per cent events of e+e~ — > J/'i> + X is from 
e+e" J/^E* -I- cc. This is a challenge for the theory, 
and also suggests the importance to measure the asso- 
ciate particles in T J/"^ + X. However, the data 
sample seems still not enough. According to the experi- 
ment d, [il], the e+e" and pairs are used to iden- 
tify and reconstruct the charmonium. Only less than 10'^ 
J/^" events are collected. Even if these events are all 
T — > J/\l/ -|- ccg, and employing the standard method 
to reconstruct all the charm pseudo-scalar and vector 



mesons that the associate charm quarks fragment to one 
can only expect to get less than 10 per cent events. On 
the other hand, it is not efficient to reconstruct all the 
important charm mesons. If only use the (— > Kir) 
channel, the number of events to be found could only be 
around 9 [3 ■ Then the statistical error may not allow a 
definite conclusion. 

However, a charm quark generally decays into strange 
quark, hence the strange particles can manifest the ever- 
production of charm. So the inclusive measurement of 
kaon may be a good probe [18]. Besides charm decay, 
the strange quark can be created from the vacuum dur- 
ing the hadronization process. Hence it is not the ab- 
solute production rate of strange particles but their en- 
hancement relative to other cases is the very signal for 
charm. For the processes we are going to investigate, 
the colour-octet models predict a flavour-singlet cluster 
(one or two colour-octet gluon(s)) associate with J/'^, 
while the colour-singlet model predicts the cluster ccg, 
which lead to strange flavour via weak decay. To find a 
measurable quantity, we study the ratio {< K > 

and < TT > are average multiplicity of and tt^, re- 
spectively) of the cluster accompanying with the J/^E". 
The cluster ccg is colour- sing let and we can straightfor- 
wardly employ JETSET [l9| to simulate its hadroniza- 
tion [2l|. We can expect the ratio is similar as those 
of e+e" — > cc — > h's at the same energy (see Table 1). 
On the other hand, it is not very practical to hadronize 
the colour-octet 2-gluon-clustcr with JETSET for the 
colour-octet process T — *■ J/vP -I- gg. Another "techni- 
cal" problem is that in the partonic level calculation, this 
cluster dominantly have a vanishing invariant mass. For 
providing a reference to the charm enhancement, we give 
the result of a cluster of gg which is in colour-singlet 
(its invariant mass is set to be the peak value of the 
ccg cluster mass). What is in common for the 2-gluon 
cluster in colour-octet and colour-singlet is that they are 
both flavour-singlet. We find that the value is 

not sensitive to the invariant mass for the colour-singlet 
gg cluster. At the same time, using the JETSET sub- 
routine (lulent), one can investigate the "independent" 
fragmentation of one gluon in the colour-octet process 
T J/^ + g. The obtained is insensitive to its 

energy. We find that the result is near that of colour- 
singlet gg cluster (Table 1). For 2-gluon system, if we 
do not care it is in colour-singlet or octet, its value 
can be obtained from independent fragmentation of each 
gluon, which is very near that of the colour-singlet 2- 
gluon cluster from the above discussion. This suggests 
that the colour-octet 2-gluon cluster should give 
not far beyond that of the colour-singlet cluster. Besides 
the above, for a comprehensive comparison, the results of 
a flavour-singlet virtual photon {M < 2mb) are also given 
as reference, but we should understand that the electro- 
magnetic interaction breaks the flavour SU(4) invariance 
(the initial quark effects). The result is also shown in 
Table 1. 

In Table 1, we give the "errors" of the calculated re- 
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e'''e~ ^ cc ^ /I's 


0.240±0.008 


99 ~* h's 


0.120±0.007 


g h's 


0.113±0.007 


7* ^ /I's 


0.176±0.008 



TABLE I: The production ratio between K and tt in various 
processes. 



suits. Their definition is based upon the following con- 
sideration: The number of collected events by the CLEO 
Collaboration is only of the order 10'^. So we should es- 
timate the statistical error, or the fluctuation of some 
10'^ event samples. Here we use rot to represent the 
"real" value of , which can be estimated by a sample 
generated by JETSET with event number large enough 
oo). At the same time, we generate many 10^ event 



samples, then we calculate 



Table 1, the values are given in the form rot ± x- 

for the processes we discuss above. We can conclude 
that the fluctuations are small and will not smear the 
difference, so that the strange particle enhancement of 
the ccg cluster is very significant, comparing with the 
flavour-singlet g (gg) as well as the virtual photon. 

In summary, we point out that the recent measure- 
ments for T decay into charmonium final states by CLEO 
Collaboration agree well with the T — > J/^* -I- ccg calcu- 
lation, provided that the value of charm quark mass is 
properly set. This indicates that the colour-singlet pro- 
cess could be the dominant for the inclusive J/'i' or ^'(25') 
production in T decay. Furthermore, by investigating the 
hadronization of the parton(s) produced accompanying 
with the J/'i' in T decay, we found that the ccg cluster 
in colour-singlet (but not flavour-singlet) will lead to sig- 
nificant strange particle enhancement comparing to the 
colour-octet (but flavour-singlet) gluon(s). Thus it is a 
good signal to probe the process T — > J/^* -I- ccg. 
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(4) 



The subscript j represents each 10'^ event sample. In 
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APPENDIX A: THE DISTRIBUTION OF ^^^^ 

<7r> 

In the following figure, we demonstrate the distribution 
of the , to show that the dispersions are really very 
small. 




FIG. 1: Here r = 



